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Determining Properties of Magnetic Elements 

^ Field of the Invention 

This invention relates to magnetic elements, particularly but not exclusive^ to 
J methods of distinguishing between magnetic elements and methods and apparatus 
for reading magnetic data tags which include one or more magnetic elements, each 
of which can differ in coercivity, saturated dipole moment Ci.e. response amplitude), 
M orientation and bias field. 

^l 

4* JO Background 

% Co-pending PCT pubUcation number WO99/35610 describes tags and reader 

^4 systems primarily intended for tags fabricated from magnetic material of low 

P coercivity, with elements at different orientations, in which data is recorded 

primarily by means of the orientation of the elements with respect to each other. 

15 The described system assumes that the coercivities of the tag elements are all the 

M same, and are very small compared to the interrogation field. 



m 



Summary of the Invention 

According to the invention, there is provided a method of reading a magnetic tag 
20 having at least one magnetic element, comprising interrogating the tag with a 

scanning magnetic field, determining transition data associated with changes in the 
magnetisation state of the at least one magnetic element, associating the transition 
data with one or more respective elements; and for each element, determining the 
element direction which corresponds to the transition data for that element. 



2S 



Preferably, the element direction is determined by selecting the direction that 
minimises the scatter of the transition field vectors resolved along the direction of 
the element. 

30 The transition data for each element can be grouped into two sets, which can be 
referred to as forward and reverse transitions. All those in the forward transition 
group have a positive component of the field vector dH/dt along the element 
vector, and all those in the reverse group have a negative component of dH/ dt 
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along the element vector. Mean field values, resolved along the element veaor, can 
be calculated. The coercivity of the element is then calculated as half the difference 
between the forward and reverse mean values, and the bias field along the element 
is calculated as the sum of the forward and reverse mean values. 

According to the invention, there is further provided a method of distinguishing 
between a plurality of magnetic elements, comprising the steps of applying a 
scanning magnetic field to the elements, determining the direction of each of the 
elements, for each of the elements, determining the components of the field in the 
direction of the element at which the element switches magnetisation states; and 
from said components, determining, for each of the elements, respective 
characteristics of the element. 

The invention further provides a method of determining, for a magnetic element, 
any one or more of a plurality of characteristics comprising the coercivity of the 
element, the local magnetic field bias resolved in the direction of the element and 
the orientation of the element, comprising the steps of applying a varying magnetic 
field to the element, determining the direction of the element, determining the 
components of the field in the direction of the element at which the element 
switches magnetisation states; and from said components, determining the one or 
more characteristics of the element. 

According to the invention, there is also provided a magnetic tag reader for reading 
a magnetic tag having at least one magnetic element, comprising means for 
interrogating the tag with a scanning magnetic field, means for determining 
transition data associated with changes in the magnetisation state of the at least one 
magnetic element, means for associating the transition data with one or more 
respective elements; and means for determining, for each element, the element 
direction which corresponds to the transition data for that element. 

Brief Description of the E>rawings 

Embodiments of the invention will now be described, by way of example, with 
reference to the accompanying drawings, in which: 
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Figure 1 is a schematic diagram of a magnetic data tag reading system; 

Figure 2 is a schematic diagram showing the components of the magnetic data tag 

reading system of Figure 1 in more detail; 

Figure 3 is a schematic diagram showing details of the signal processor/ controller 
5 illxistrated in Figures 1 and 2; 

Figures 4 and 5 illustrate the receive coil set; 

Figures 6 and 7 illustrate the transmit coil set; 
M Figure 8 illustrates the antenna comprising transmit and receive coil sets; 

Figure 9 is a flow chart illustrating the overall processing algorithm; 
41 10 Figure 10 is a schematic diagram of data acquisition circuitry; 

1^ Figure 11 illustrates the transmit current waveforms; 

[ij Figure 12 illustrates a flow chart for the signal processing and filtering algorithm; 

^ Figure 13 illustrates signals at the inputs to the ADC from the x, y and z receiver 

B 

qj coil preamplifiers, for a single element transition; 

^1 15 Figure 14 illustrates the composite filter output of the signal in Figure 13; 

'y, Figure 15 illustrates a 3D scatter plot of the filtered receiver vectors; 

P Figure 16 illustrates a flowchart for the clustering algorithm used for planar tags; 

* ' Figure 17 illustrates the composite filter output for three parallel elements; 

Figure 18 illustrates a flowchart for a parallel element clustering algorithm; 
20 Figure 19 illustrates a 3D scatter plot of the transition field vectors for a single 
element; 

Figure 20 illustrates the same 3D scatter plot as Figure 19, tilted such that the 
transition planes are edge-on; 

Figure 21 illustrates distribution of field vectors that occur along a misaligned 
25 element direction vector; and 

Figure 22 illustrates a flowchart for the calculation of the mean switching field, 
switching field variance, coercivity and DC bias field. 

Detailed Description 
30 Referring to Figure 1, a magnetic tag reading system comprises a magnetic data tag 
' 1, an interrogation unit 2 and a signal processor/controller 3. Magnetic tags 1 to be 
used with a magnetic tag reader according to the invention can record information 
by means of elements of differing coercivities, local bias fields and response 



wo 01A3321 



-4- 



PCT/GBOO/03092 



amplitudes, as well as orientation. This includes tags described in PCT publication 
number WO99/35610, as well as tags described in, for example, US 5,204,526, US 
5,729,201 and W098/26312. In general terms, magnetic tags 1 comprise magnetic 
elements which typically switch magnetisation state, for example magnetisation 
direction, at given values of applied field depending on element properties, for 
example coerdviiy. These elements include, for example, thin film elements, 
bistable elements, Barkhausen wire elements and high-permeability elements. The 
applied field which causes switching depends on the magnitude of the component 
of the interrogation field vector in the direction of the element. 

Referring to Figure 2, the tag 1 is attached to an item being labelled or tagged 4, and 
is placed within an interrogation volmne 5 within the interrogation imit 2. The 
interrogation unit 2 includes an antenna 6, which comprises transmit and receive 
coil sets 7, 8. The tag 1 is interrogated by a scanning magnetic field 9 generated by 
the transmit coil set 7 under the control 10 of the processor/ controller 3. In 
response to the interrogating magnetic field 9, the tag 1 generates a detectable 
magnetic field response 11, which is detected by the receive coil set 8. The 
processor/controller 3 receives input signals 12, 13 from the transmit and receive 
coil sets 7, 8 respectively and processes the signals to decode data stored on the tag, 
which is made available at an output 14. 

Referring to Figure 3, the processor/controller 3 comprises a waveform generator 
15 for driving the transmit coil set 7, data acquisition circuitry 16 for receiving 
respective input signals 12, 13 from the transmit and receive coil sets 7, 8 and a 
digital signal processor 17 for processing the resulting output signals 18 from the 
data acquisition circuitry 16 to provide the decoded tag data 14. 

The transmit and receive coil set arrangement 7, 8 is described in detail by reference 
to Figures 4 to 8. 

Figures 4 and 5 illustrate the receive coil set 8. The receiver coils are constructed 
on a cylindrical former 20 of diameter 200 mm and length 400nun. Figure 4 
illustrates the three sets of orthogonal coils used to couple with the tag magnetic 
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elements within the interrogation zone. For the y-direction, the receiver coil set 
comprises 4 coils 21, 22, 23, 24. Inner coik 22, 23 lie on the former 20 and extend 
120mm along the x-direction 25. Both inner coils 22, 23 comprise 100 turns 0.4mm 
ecw. The outer coils 21, 24 comprise 58 turns of 0.4mm ecw and are woimd on a 
second co-axial former (not shown) 260mm in diameter. The coils extend 156mm 
along the x-direction. The four coils 21, 22, 23. 24 are connected in series in the 
electrical sense illustrated and 'balanced* by small mechanical re-alignments to 
achieve zero sensitivity to a uniform magnetic field. A second receiver coil set as 
illustrated is sensidve to tag generated field in the z-direction. This coil set is 
identical to the coils 21, 22, 23, 24 but rotated through 90^ as shown. The third coil 
set sensitive to tag generated field in the x-direction comprises two solenoid coils 
26, 27. The inner coil 26 comprises 100 turns 0.4mm ecw wound on the former 20, 
and is 120mm long. The outer coil 27 comprises 58 turns of 0.4nun ecw wound on 
the second 260mm diameter co-axial former and is 156mm long. Figure 5 illustrates 
all the coils wound on the inner former 20, and the outer former 28. 

Figure 6 illustrates the three orthogonal transmit coils configuration 7. The coils 
are wound on a cylindrical former 30, 370mm long and 300mm diameter. A 
uniform magnetic field in the y-direction is produced by four coUs 31, 32, 33, 34. 
First and third coils 31, 33 comprise a 'modified Helmholz' arrangement similar to 
coils 15 and 16. Second and fourth coils 32, 34 comprise a second modified 
"Hehnholz' arrangement, with a magnetic axis 25° offset from the first and third 
coils 31, 33. The two 'modified Hehnholz' coil sets have magnetic axes 12.5° either 
side of the y-direction. The first coil 31 comprises 50 turns 1.4mm ecw and extends 
370mm in length along the former. Where this coil 31 connects across the open 
end of the former 30, the coil is a flaaened half circle with the total coil aperture 
width of 570mm. The two edges of the coil 31 that lie along the solenoid (x- 
direction) subtend 120° at the axial centre of the former. Second to fourth coils 32, 
33, 34 are identical in size and form. Their orientation around the former 30 is 
described above. The four coils are connected in series in the sense illustrated. A 
second transmit coil set generates uniform field in the z-direaion. This set 
comprises four identical coils orientated in an orthogonal direction as illustrated. 



wo 01/13321 PCT/GB00/O3OW 



The final transmitter coil consists of a long solenoid coil 35 comprising 260 turns of 
1.4mm ccw on the coil former. This generates uniform field in the x-direction. 

Figure 7 shows the overall transmit coil arrangement for generating imiform field in 
5 three orthogonal directions. Figure 8 illustrates the antenna 6. The transmit coils 
on the former 30 are located co-axially with the receiver coil tube 20. The 
interrogation volimie 5 is defined by a further 190mm ID co-axial tube (not shown) 
that is used to define a mechanical constraint on possible tag positioning in the 
antenna 6. The longitudinal region of highest sensitivity is less than 10cm long and 
10 tags can be accurately read when separated by 10cm or more along the axis of the 
reader tube. 

Figure 9 illustrates the overall sequence of steps required to decode data stored on a 
magnetic tag. The first stage is data acquisition (step si). Data is acquired by 

25 detecting the field 1 1 resulting from the application of a scanning interrogation field 
9 to the tag 1, digitising the resulting signals and storing them for subsequent 
processing. This results in 3 channels of input data, one for each of the x, y and z 
directions. Digital signal processing is carried out to identify individual switching 
points, also referred to herein as transitions (step s2). This results in an array of 

20 transition information. Each transition is associated with an element (step s3) to 
provide an array of elements. The elements are then individually decoded (step s4). 
Finally, the tag is decoded to provide tag value data (step s5). 

Data Acquisition 

23 Referring to Figure 10, an example of the signal processor/controller 3 according to 
the invention comprises a National Instruments PCI6711 4-channel DAC card 39 
for waveform generation and a National Instruments PCI6110E 4-channel ADC 
card 40 for data acquisition. The cards are mounted into an industry standard IBM 
compatible PC 41 running Windows 95™. The waveform generation card 39, under 

30 software control, generates three transmit excitation voltages 42, 43, 44 which are 
passed through respective low-pass filters 45, 46, 47 and amplified by respective 
power amplifiers 48, 49, 50 to drive respective orthogonal transmit coils, which are 
arranged in a series resonant configuration with respective capacitors 51, 52, 53 and 
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resistors 54, 55, 56. The drive current is for example 3A rms to generate a 2.5kA/m 
interrogation field. The transmitter currents are monitored by respective current 
sense resistors 57, 58, 59 which are fed through respective amplifiers 60, 61, 62 as 
inputs 63, 64, 65 to the data acquisition card 40, where they are digitised at a sample 
rate of, for example, 160kH2. The instantaneous transmit field vector can be 
determined from these three signals with knowledge of the relationship between the 
transmit coil field and the current response. For example, pre-calibration of the 
system is carried out by measuring the transmit field for different values of driving 
current. 

Signals induced in the orthogonal receive coils are amplified by respective amplifiers 
66, 67, 68, filtered by respective 130H2 notch filters 69, 70, 71 to remove any 
transmit field component, and fed as inputs 72, 73, 74 to the data acquisition card, 
where they are digitised at a sample rate of, for example, 160kHz. The data 
acquisition is buffered in such a fashion that data is clocked into a buffer, and read 
from the buffer asynchronously at some later point. The buffer depth is sufficient 
to accommodate the worse-case latency in the subsequent processing step. 

A continuous scan is used in this example to interrogate the tag, based on a nominsJ 
130 Hz rotating magnetic field, whose normal vector is arranged to trace out a spiral 
scan over the surface of a complete sphere, tracing a path from one pole of the 
sphere to the other and back. The equations for the components of the 
^transmitted* B interrogation field are given by: 

« (cos^(<l)) *cos(e) + sin^ (<{>)) *cos (cor) 
+ (sin(i|))*cos((t>)*cos(e) - cos((|))*sin((|)))*sin((or) 

By - (cos(<j))*sin(<|>)*cos(e) - sin(<|))*cos(ij>))*cos(cor) 
+ (sin^ (<j))*cos(e) + cos' (*))*sin(a)t) 

B^ - (-cos((|>))*sin(e)*cos(cDr) - sin(e)*sin(<|>)*sin(©t) 
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where t is the time, © is the angular frequency of the 130Hz scan, ^ - (constant)* 9, 
and 9 » cos**(l-t/7). 7 is the total time for one complete interrogation. 

Figure 11 shows the three transmit current waveforms 80, 81, 82 received at the 
5 inputs 63, 64, 65 of the data acquisition card 40. 

Digital Signal Processing 

The digital signal processing stage (step s2) performed on the data input to the data 
acquisition card 40 is now described with reference to the flow chart description of 
10 the processing algorithms in Figure 12. 

The purpose of the DSP algorithm is to identify individual transitions, and to record 
all the relevant parameters on each transition for subsequent processing algorithms. 
This leads to a large reduction in the volume of data passed on to the subsequent 
IS processing stages. 

The DSP algorithm operates on the three channels of sample data produced by the 
data acquisition process. Figure 13 shows the raw impulse responses in the x, y and 
z channels for a single transition. 

20 

Referring to Figure 12, in a first step slO, an FIR filter is applied to all three 
channels to produce three sets of filtered data, which form a receiver vector. The 
simplest filter consists of three rectangular sections, and provides a method of 
measuring the height of the peaks in the raw data. If the central section has width 
2S w and height + 1, then the outer two sections have width w/2 and height -1. The 
width, w, is typically the same value as the response time of the magnetic element, 
for example, 20-30p,s. 

The transmit field vector, H, is used to determine the correct polarity of the element 
30 transition in each of the x, y and z receiver coils. The transition polarity in any 
receiver coil is directly related to the polarity of the rate of change of the field 
vector dH/dt in the direction of the receiver coil. dH/dt values are used to 
produce a "polarity vector", where each component can take the value ±1. The 
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scalar (dot) produa of the polarity vecxor with the filtered receiver vector is 
calculated (step sll) and the receiver vector magnitiide, a positive number, is 
multiplied by the sign (±1) of the result (step sll). This results in the composite 
signal shown in Figure 14, in which the polarity of transitions for every element is 
always the same, allowing the use of a simple peak detector to determine peak 
values. 

Peak detection techniques are well known in the art. In this case, a simple threshold 
is used to gate the peak detector input data, to avoid noise appearing as spurious 
peaks. A peak is identified when three or more values exceed the threshold (step 
sl3), and where the current value is greater than both the previous and next value 
(step sl4). The time of the peak is interpolated to a greater resolution than the 
sample frequency by a simple quadratic fit to these three points. 

The data for each transition is stored in an array (step si 5). The data includes: 
Time 

Field vector (H) 

Rate of change of field vector (dH/ dt) 
Receiver vectors (both raw and FIR filtered) 

Element Association 

The function of the element association algorithm is to associate transition data 
points with particular magnetic elements in the tag. Subsequent processing steps 
can then analyse the data for each magnetic element in isolation, thereby reducing 
an apparently complex problem with multiple elements into a series of relatively 
simple numerical solves. 

There are two primary mechanisms that are used to associate transitions with 
elements, depending on whether the elements are generally parallel or not. These 
are described below. In the general case, the first step is to separate into groups 
using a non-parallel algorithm, and then, if required, to analyse each separate group 
to see if it contains more than one parallel or near-parallel elements. 
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For non-parallel elements, the filtered receiver vectors are used to separate out the 
transitions between elements. This can be clearly seen from Figure 15, which shows 
the filtered receiver vectors in a 3D scatter plot for an example tag having seven 
non-parallel elements. Inspection of the plot shows that the majority of the 
transition points lie along one of 7 different lines through the origin, which 
indicates that there are 7 discernible directions of elements in the example tag. 
Each direction can be described by two parameters, and therefore the transitions 
can be clustered together into groups in 2D. There are a nxmiber of different 
appropriate techniques than can be used to achieve this multi-dimensional clustering 
(e.g. S. Makeig, S. Enghoff, T-P. Jung, M. Westerfield, J. Townsend, E. Courchesne 
and T.J. Sejnowski, "Moving-Window Independent Component Analysis of Event- 
Related EEG Data: Component StabiUty, Journal of Neurophysiology''). Additional 
knowledge about the particular tag construction can be useful to simplify the 
problem. For example, if all the elements are in the same plane, then the problem 
can be reduced to a one-dimensional problem. Knowledge of the number of 
elements expected can assist in the clustering process. 

In the particular case of a planar tag, with a known number of elements, the 
algorithm outlined in Figure 16 is used. The normal to the plane of the transitions 
is determined by, for example, a numerical process (step s20). For example, the dot 
product of every receiver vector with an estimated direction vector is calculated, 
and this process is iterated until the sum of the magnitude of the dot products is 
minimised. This reduces the problem to a 1-D problem i.e. the angle in the plane. 
An in-plane set of vectors can be calculated from the original set of vectors simply 
by subtracting from each vector in turn the dot product of itself with the normal to 
the plane. In-plane angles between any two in-plane vectors can them simply be 
calculated in the usual way using dot products. All the in-plane angles are wrapped 
into the range 0-180** by adding or subtracting mxiltiples of 180** as required. The 
algorithm calculates a histogram of in-plane angles relative to some arbitrary datima, 
such as the first point. For example, if the histogram bins are 1** wide, then the nth 
bin will contain a count of the number of angles that fall in the range to (n+ 1)**. 
This will typically give a series of peaks, one for each element. For example, the 
algorithm obtains the second point from the transition array (step s21), measures 
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the in-plane angle relative to the first point (step s22) and increments the 
appropriate bin of the histogram (step s23). This process is repeated until all the 
data has been processed (step s24). After applying Gaussian smoothing to the 
histogram data (step s25), the direction of an element in the tag can be fo\md by 
determining the highest peak in this histogram (step s26). To determine the 
transitions that belong to the element in a given direction, the algorithm finds all 
the transitions that are within, say, 2** (m plane) of this direction (step s27). The 
processed peak is then eliminated from the calculation (step s28) and this 
processing sequence is repeated (steps s26 to s28) imtil all the data has been 
processed (step s29). 

To separate parallel elements, the algorithm makes use of two properties of a 
continuous scan of the field vector, H, aroxmd the elements, first that the elements 
transition in order from the lowest to the highest coercivity and second, that the 
field vector, H, rotates by at least 90'' between the last transition of the highest 
coercivity element in one direction and the first transition of the lowest coervicity 
element in the reverse direction. 

If some of the elements do not change state, because the transmit field does not 
reach a high enough value, then there will be fewer transitions, in a 180*" scan, than 
there are elements. In this case, transitions are "lost", starting with the highest- 
coerdvity element. Figure 17 shows the filtered composite waveform over a few 
rotations of transmit field, for three parallel elements with different coercivities. 
Each element is associated with a respective peak 90, 91, 92 and it is relatively 
straightforward to separate out the transitions belonging to different elements. An 
outline algorithm to achieve this is shown in Figure 18. 

This works by maintaining an element counter that is incremented each time a new 
transition is identified, and set to zero each time the field rotates by more than 90** 
between transitions. Data is extracted from the transition array (step s30) and the 
algorithm determines whether the transmit field has rotated by more than 90** since 
the last transition point (step s31). If it has, the element coimter is reset to zero 
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(step s32)- Following this, the element niunber corresponding to the transition is 
set according to the current value of the element counter (step s33). The element 
counter is then incremented (step s34) and the process repeated for the next 
transition point (step s30). For example, the first transition following the zeroing of 
the element coimter is associated with element 0, the next with element 1 and so on, 
until the field has rotated by more than 90 degrees. 

Element Decode 

The purpose of the element decoding algorithm is to take transition data belonging 
to one element, and to determine the best-fit direction vector for this element. 
Once the direction is known, the coerdvity of the element, and any net DC field or 
'"bias'' along the element vector can be calculated. 

Figure 19 is a 3D scatter plot of transition point field vectors, for a single bistable 
magnetic element with a finite coercivity. In this example, the field has been 
scanned approximately over the surface of a sphere, so the transition points lie 
roughly on two circles 93, 94. More generally, the transition points would be 
expected to lie on one of two planes. By lilting the view of the transition data in 
the scatter plot, it is possible to show the two planes edge-on, as illustrated in 
Figure 20. The bold vertical arrow 95 shows the element vector. 

The element decoding algorithm attempts to determine the best vector direction for 
the element, by minimising the scatter of field vectors resolved in this direction. 
Figure 21 illustrates the situation where a guess 96 has been taken for the element 
veaor that is not in the correct direction. Taking the upper set of transitions 93, it 
is clear that when these are projected onto the element vector 96, they form an 
extended distribution 97 (shown by a darkened section) along the vector 96. As the 
vector is rotated aro\md, the extent of this distribution will be smallest when the 
vector is closest to the actual direction of the tag element. 

Figure 22 shows a flowchart for the algorithm used to calculate the error used in the 
iterative solving process. The algorithm uses the current guess for the element 
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vector direction, V. Initially, this is the vector direction from the element 
association algorithm. 

The first data point is retrieved (step s40) and the dot product of the direction 
estimate V with the polarity vector for dH/dt is calcxilated, as described in relation 
to Figure 12 above (step s41). If this is positive (step s42) i.e. for the upper set of 
transitions, then the dot product of the field vector and the element vector is 
calculated (step s43) and added to the upper set of statistics (step s44). The dot 
product resolves the component of the field vector along the element vector. The 
same calculation is carried out for the lower set of transitions (steps s45, s46), 
indicated by the negative dot product at step s41. The upper and lower (forward 
and reverse) sets of transitions are distinguished by the sign of dH/dt along the 
direction of the element, or alternatively by the sign of the fikered receiver vector 
along the direaion of the element. This procedure is repeated for all the data 
points (step s47). An average value of variance is calculated from the variances for 
each of the upper and lower sets of transitions separately, weighted by the number 
of transitions in the upper and lower sets of transitions (step s48). A standard 
formula is used to calculate the variance for each set of data. For a set of measured 
data points, x, the variance, var(x) is the mean of the squares of x minus the square 
of the mean of x, or mathematically 

var(x) = (x^)-(x)^ 

The weighted variance of //„ upper transition points, «, with variance var(ii) and 
Ng lower transition points, /, with variance var(/) is then given by: 

var(w) + //, var(/) 
var = • 



The weighted variance is used as a measure of the error in the guessed vector 
direction. When the guessed direction is equal to the actual element direction, the 
weighted variance will generally have its minimum value. The value will never fall 
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to zero, because there is always a certain amount of noise in the determination of 
the transition field, arising from sources such as electronic noise and randomness in 
the material behaviour. In the simplest case, the variance is a function (nimierically 
evaluated, rather than an analytic function) of two direction variables, such as 9 and 
4» from the spherical polar co-ordinates (r, 9, <}>). The value of this function can be 
minimised using a standard numerical minimisation algorithm. The variance varies 
approximately quadraticalfy with the deviation from the ideal direction, and this 
means that the minimisation algorithm can be extremely efficient (the "quadratic" 
case is generally considered to be the easiest). Multi-variate numerical minimisation 
algorithms are well known in the art - for example, Powell's method, 

When the weighted variance is not minimum (step s49), the direction estimate V is 
adjusted in accordance with the appropriate minimisation algorithm (step s50) and 
the algorithm re-run with the new value of V. When the weighted variance is 
minimised, the mean values of the field for the upper and lower sets of transitions 
are calculated (steps s51, s52). The coercivity of the element is calculated as half 
the difference between the two switching fields (step s53), while the DC field along 
the element is calculated as the sum of the two switching fields (step s54). 

Additional parameters as well as direction can usefully be added to the nimierical 
minimisation. The most important term to add is the vector velocity, which allows 
the algorithm to deal with the movement of the tag elements during the decoding 
process. The element direction in steps s41, s43, s45 is then a function of the time 
at which the transition occurs, and the funaion for the variance then depends on 
four parameters (for example 9, <t>, d9/dt and d<|>/dt). Once again, this function can 
simply be minimised using a standard multi-variate minimisation algorithm. 

Many magnetic elements do not behave ideally, and show significant changes in 
their switching field (or coercivity) depending on the value of dH/ dt. A general 
form for the switching field, H^i^, resolved along the element is: 
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where a,b etc arc arbitrary powers. If the coefficients k are known, then the value 
of Ho can be calculated from the measured switching field, and the variance of this 
value can be minimised as before. If the coefficients k are not known, but the 
values a, b are known, then the function for the numerical minimisation can also 
5 include the coefficients, k, as function arguments, as well as the direction and 
velocity terms. In this case, the values of k can be used to distinguish between 
different types of materials and thereby store more data. 

Anisotropic thin-film magnetic materials can exhibit a further form of non-ideal 
10 behaviour. For materials with an easy-axis of magnetisation, the in-plane field 

perpendicular to the easy axis can influence the field at which the material switches. 
A similar approach to the one described above can be used to calculate a nominally 
constant value, Ho, from the raw switching points. 

15 Tag Decode 

The primary output data for each magnetic element from a reader according to the 
invention is as follows: 

Orientation in the reader (vector) 
20 Coercivity of the element (scalsir) 

Bias field along each element (scalar) 
Amplitude response (scalar) 

Secondary data for each element includes 

dH/dt coefficients 
Perpendicular field coefficients 
Response time 

Characteristic response "shape** or spectrum 
30 Statistical distribution of primary parameters 

This data assumes little about the construction of the tag. The structure of the tag 
(e.g. which elements share bias magnet elements) may be used to provide more 
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detail - for example, the magnitude and direction of an overall bias field. The 
details of the chosen coding scheme are used to translate all these raw parameters 
into useful data stored on the tag. 

The above examples of the invention are intended to be illustrative, rather than 
restrictive. A person skilled in the art would imderstand that various modifications 
and variations in the detailed implementation are possible, and are considered to be 
within the scope and spirit of the invention as defined in the appended claims. 



